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In  Ref.  1  Chapman  considered  the  problem  of  laminar  mixing  at 
constant  pressure  for  a  fluid  with  Prandtl  number  one,  and  a  viscosity 
power  law  of  the  form  ^  T^.  Because  of  the  recent  interest  in  the 
flu id -mechanic  description  of  hypersonic  wakes,  mixing  problems  are 
being  again  investigated  intensely  both  experimentally  and  theoreti¬ 
cally.  As  an  example.  Ref.  2  extends  Ref.  1  by  assuming  a  Blasius 
starting  velocity  profile  rather  than  a  uniform  one.  (A  uniform  pro¬ 
file  was  assumed  by  Chapman  as  a  necessity  imposed  for  the  conserva¬ 
tion  of  similarity) . 

It  will  be  shown  in  this  brief  note  that  the  velocity  along  the 
dividing  streamline  can  be  determined  in  a  simple  manner  by  approxima¬ 
ting  the  Integral  solution  of  Ref.  1.  The  effect  of  finite  base  rad¬ 
ius  (or  Reynolds  number)  is  also  investigated  using  this  approximate 
solution. 

In  terms  of  the  stream  function  the  differential  equation  of 
motion  is: 

2  dC  dC  V®  dC  ° 

where 

C  =  Y/VUvsC,  g(C)  =  uT^"^ 


T  and  u  are  nondimen sional  with  respect  to  the  free-stream,  the  sub¬ 
scripts  d  and  o  refer  to  the  "dead  water"  region  and  stagnation  point. 
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and  M  is  the  free-stream  Mach  number.  Chapman's  boundary  conditions 
are: 

At£=oo:  u=l,  C  =  u  =  0  (2) 

The  velocity  at  the  dividing  streamline,  Uj^,  is  given  after  a 
formal  solution  of  Eq.  (1)  as  follows: 


I  m 


I  FdC  +  J  FdC 


•CO  O 

where 

«p  ^  'ic} 

F«)  -  - 1 -  (4) 

We  make  the  observation  that  in  the  expression  for  the  Integrand 
F,  the  contribution  of  the  exponential  part  in  the  numerator  is  stronger 
than  the  denominator  for  Increasing  g  (or  Q) .  Hence  in  the  intervals  1: 
-00  <  2  ^  0  and  II;  0  ^  ^  +  oo  most  of  the  contribution  comes  about 

from  the  values  of  g  corresponding  to  the  highest  Q  inside  the  inter¬ 
val.  For  a  first  iteration  let  us  therefore  assume  for  u  the  following 
step  function:  Inside  I,  u  =  u^;  Inside  II,  u  >  1.  Simple  Integration 
of  Eq.  (3)  yields 


Assuming  that  and  u)  >  0.75  the  calculations  show  that  for  M  *  0, 

1,  2,  3,  4,  5,  7,  10,  15,  20,  the  corresponding  values  of  u^^  are: 

0.570,  0.573,  0.581,  0.590,  0.598,  0.605,  0.619,  0.634,  0.652,  0.664. 

(3) 

For  M  >  0  and  5  Chapman'  '  gives,  through  an  exact  numerical  solution, 
u^  -  0.587  and  0.597.  Comparison  shows  that  our  closed  form  approxima¬ 
tion  Is  In  error  of  less  than  -3%  and  +1.5%  correspondingly. 
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(4) 

Experiments  show  '  that  Uj^  is  a  function  of  Reynolds  number.  The 
analysis  of  Ref.  2,  in  which  the  influence  of  an  initial  finite  bound- 
ary  layer  thickness  was  studied,  through  the  parameter  s  ~ 
yields  results  which  are  independent  of  the  base  radius  r^.  (This  oc¬ 
curs  because  as  it  can  be  seen  from  Fig.  1,  r  =8,  sin  a  =  s  sin  P) . 

o  D  n 

One  might  conjecture  that  one  way  to  introduce  the  Reynolds  number 
would  be  to  assume  that  u  =  0  not  at  £  =  -oo  but  at  £  =  where  is 
finite  and  positive.  This  assumption  implies  a  finite  radius  in  the 
direction  perpendicular  to  the  main  flow.  Following  the  same  method 
used  for  the  derivation  of  Eq.  (5)  we  find; 


Erf.  (X^) 
Erf.  (X„)  + 


(6) 


where  and  Erf.  denotes  the  error  function.  Figure  2 

shows  the  function  for  three  different  Mach  numbers. 

These  results  are  best  interpreted  in  the  physical  plane  s,y. 
For  simplicity  assume  M  =  0  so  that^^^ 


(7) 


Let  correspond  to  y^  and  s^  where  the  subscript  n  indicates  the 
position  of  the  "neck"  of  thickness  h  as  shown  in  Fig.  1.  To  calcu¬ 
late  the  Integral  in  Eq.  (7)  we  approximate  u(Q  in  the  interval 
0  ^  ^  by  dropping  the  first  term  in  Eq.  (1) .  A  simple  integra¬ 

tion  yields: 


u 


(8) 


Introducing  the  above  into  Eq.  (7)  we  have: 
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Fig.  1  ~  Schematic  diagram  of  base  flow 


Fig.  2 —The  veixity  at  the  dividing  streamline  for  different 
boundary  conditions  inside  the  base  region 
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Assuming  that  all  of  the  mass  contained  in  the  boundary  layer  over  the 
body  goes  through  the  neck,  we  set 


The  f'ctor  40  is  fixed  by  the  experiments  of  Ref.  3.  Setting  the  posi¬ 
tion  of  the  neck  \  radii  back  of  the  base,  Eq.  (9)  yields,  after  some 
trigonometry: 


=  20  u. 


'k  cos 


From  this  equation  it  appears  that  the  value  of  Uj^  is  again  directly 

Independent  of  the  Reynolds  ntunber.  Since  the  angle  P  is  of  the  order 

of  10°,  so  that  cos  P  »  1,  and  since  from  experiments,  \  is  of  order 

one,  the  assumption  C  «  is  still  reasonable.  The  proof  of  this  last 

n 

statement  lies  in  the  fact  that  Eq.  (11)  yields  a  value  of  =  12 
when  Up  »  0.6,  X,  w  1  and  cos  g  ~  1;  then.  Fig.  2  indicates  that  our 
choice  for  Up  is  consistent  with  the  fact  that  at  ■  12  the  asymp¬ 
totic  value  for  Up  has  been  reached.  In  fact  it  is  reached  roughly 
when  Cjj  ^  3.0.  One  needs  a  neck  length  of  the  order  of  10  radii  in 
order  to  make  a  correction  in  Up  for  finite  base  radius.  However, 
such  values  are  not  observed  experimentally. 
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